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Abstract 

We propose novel types of parametric oscillators generating both three- 
photon and four-photon bright light which are accessible for an experiment. 
The devices are based on the cascaded down-conversion processes and consist 
of second-order media inserted in two-resonant mode cavity. Discussion of dis- 
sipation and quantum features of the system are performed by the quantum- 
jump simulation method and concerns to the Wigner functions. The phase- 
space multistabilities and critical threshold behavior of three- and four-photon 
subharmonics are obtained. 

PACS number(s): 42.50.Lc, 42.65. S, 42.50.Dv, 42.65.Yj 
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I. THE PROBLEM 



Optical parametric oscillators (OPO's) based on processes of down- conversion in a cavity 
have proved to be efficient sources of light with a range of unique properties. Among them 
it should be noted the nonclassical nature of light, including quantum correlations, sub- 
Poissonian statistics and a larger amount of squeezing JD. OPO provides a method for 
producing tunable radiation from a coherent fixed-frequency pump source in the intracavity 
down-conversion, where a pump-photon splits into a pair of subharmonic photons. So far 
OPO has only been realized experimentally for two-photon down-conversion in nonlinear 
media with second-order x susceptibility. Experimental and theoretical results can be 
found in [pf. 

In this letter we propose a new conception for the realization of multiphoton parametric 
oscillators using the idea of composite multiphoton interaction. Two schemes of OPO based 
accordingly on a three-photon and four-photon intracavity down-conversion are presented. 

One of the principal motivations for the investigation of multiphoton OPO concerns to a 
production of nonclassical states of light with the novel features. Other motivation is their 
possible exploitation for quantum information technologies, because sources offering a great 
variety of multiparticle entangled states, first of all meaning photons, are required for the 
implementation of many quantum communication processes such as direct transmission of 
the state, quantum teleportation, quantum cryptography, and entanglement swapping ||. 

It is known that two-photon down- conversion so far is the most standard source of pair 
of polarization-entangled photons which is relevant to many problems of quantum informa- 
tion (see, previous experiments on photon splitting in ||, and for polarization-entangled 
photon pairs ||). Although two-particle entanglement has long been demonstrated experi- 
mentally the realization of entanglement between three or more particles has been consid- 
ered as a difficult problem. Among the possible sources of three-photon states (so-called 
Greenberger-Horne-Zeilinger states ||) we mention the process of degenerate three-photon 
down-conversion co^u/3-\-uj/3 + uj/3 in a medium. This process has been studied for 
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running waves 0. Note, that most recently the polarization entanglement for three spatially 
separated photons has been observed ||. Three-photon down-conversion in an optical cavity 
has also been proposed and investigated in |J , but too weak interaction between photons in 
•^(3) medium makes experimental realization of this scheme a difficult problem. As regards 
to a four-photon splitting uj-^lu/A + uj/A + uj/A + uj/A in a x medium, as we know, this 
process was not considered yet. 

This letter reports the analysis of two schemes of multiphoton OPO in a double resonant- 
mode cavity involving two crystals of second-order susceptibilities. They are the following: 

A. Four-photon splitting in cascaded down-conversion 

The scheme proposed contains two sequential nonlinear processes of degenerate two- 
photon splitting in x^ media taking place inside the same cavity. This scheme for the 
special configuration when the coupling in- and out-fields occurs at one of the ring cavity 
mirrors is shown in Fig.l. The pump frequency uj is converted to the subharmonic frequency 
co>2 = uj/2 and then into the frequency uj\ = uj/A in the cascaded parametric processes 
uj — > u>2 ~\~ u>2 and uj 2 — * Wi + ^i ■ We assume that the each of these processes takes place 
effectively only in the definite nonlinear crystal. The cavity is also assumed transparent at 
the pump frequency. So, we neglect the pump depletion effects considering the amplitude 
of driving field as a classical constant. We model the combination of these processes by the 
interaction Hamiltonian: 

Hi = il%Xi{Ea} 2 2 — E*a\) + ihki(a^a 2 — a^a^), (1) 

where a± and a 2 are the operators of the modes uj\ and 102-, Xi an d k\ are the nonlinear 
coupling constants for the processes uj —>■ 102 + ^2 arid LO2 — ^ <-^i ~t~ 001 respectively which 
are related to second-order nonlinear susceptibilities x^ ■ E is proportional to the coherent 
driving field amplitude. 
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B. Three-photon splitting in cascaded down-conversion 



This system mainly is similar to the previous one. The only difference is that the pump 
frequency u converts to two different frequencies u>i = uj/3 and u 2 = 2u/3 in the parametric 
interaction uo = io\ + 002- The subharmonic u)\ is also produced in the other process of down 
conversion 002 = u)\ + U\. On the whole two subharmonic modes U\ = u>/3 and uj 2 = 2u;/3 
are created in the cavity due to the cascaded processes. In the undepleted pump field 
configuration this model is described by the following interaction Hamiltonian: 

H\ = ihx2{Ea\a,2 — E*a\a2) + ihk2(a\ 2 a2 — o^ai,), (2) 

where a± and ai are the operators of oj\ and UJ2 subharmonics, and \i an d ^2 are the coupling 
constants of the processes uj — » u)\ + u 2 and to^ — > W\ + uji, respectively. 

We note, that optical schemes involving cascaded nonlinearities in a cavity are currently 
of considerable interest and hold promise for applications including, in particular, frequency 
tunable sources of light, large third-order nonlinear effects via cascaded second-order non- 
linearities and squeezing phenomena. The references and experimental observation can be 
found in jnj. Possible operational regimes of a combined system, as a rule, essentially differ 
from those in pure processes in each of nonlinear media. Some results in this direction are 
devoted to the cascaded frequency doubler [IT]. Concerning our schemes, it is interesting 



to note, that the stable stationary above-threshold regime of oscillation is realized even in 
undepleted pump regime. The other novelty here is comparatively low values of the thresh- 
old of generation due to the schemes proposed can be accessible for experiments. Contrary 
to polarization-entangled state generation in two- and three-photon devices [§,|J the basic 
multiphoton states of our systems are neither based on polarization nor on momentum, but 
on phase-space entanglement of modes. 
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II. FOUR-PHOTON OPO 



The system of interest is dissipative because the subharmonic modes suffer losses due to 
partially transmission of light through the mirrors of the cavity. So that the reduced density 
operator p of both modes obeys a master equation 
dp 1 

— = — [ifi, p] + ^ li( 2a iP a l - a UiP ~ pa\ai), (3) 

i=l,2 

where 71 and 72 are the cavity damping rates for the modes lo\ and uj 2 . 

To proceed with further analysis we present now the semiclassical steady-state solutions 
and stability properties of the system proposed. In the classical limit the system may be 
described by the complex-field amplitudes a\ and a 2 of the resonant modes lo\ and uj 2 whose 
real and imaginary parts respectively represent the dimensionless position and momentum. 
They obey the equation of motion: 

— — = -7^1 + lk X OL 2 OL^ 

at 
dot 

— = -72« 2 + 2xi«2^ - haj. (4) 

Solving these equations for the steady-state regime and carrying out the standard linearized 
stability analysis one can arrive at the following results. 

The trivial zero- amplitude solution a.\ = a 2 = is stable in the region E < E th , 
Eth = 72/2x1, with E t h being the threshold value of the pump field amplitude, and it 
describes the below-threshold regime of oscillation of both subharmonics. The stable above- 
threshold solutions, for E > E t h, expressed in terms of photon numbers rij and phases ipj of 
subharmonic modes (ctj = exp(iifij)y/rTj, j = 1,2) are calculated as: 

m = lll2 (E/E th - l)/2k\, n 2 = -yf/Akl, (5) 

$ 7T $ 

ifi = — ± -m, ip 2 = — ±im, n,m = 0, 1,2..., (6) 

where $ is the phase of the driving field. 

So, this scheme of multiphoton OPO allows us to reach stable generation of subharmonics 
at frequencies lo\ = cu/A and uj 2 = uj/2 at E > E th in undepleted pump-field regime. It 



should be pointed out that the occurrence of such a regime of oscillation becomes possible 
just due to the combination of two kinds of down-conversion processes in a cavity. One of 
the peculiarities of cascading dynamics is that the threshold pump power P t h = chwE^ h /2L 
(L is the optical way of pump field) is only related to the coupling constant xi °f u 
u/2 + u/2 parametric process. This coupling constant is proportional to the second-order 
susceptibility that leads to the comparatively low value of the threshold. We illustrate this 
fact for the system depicted in Fig.l, choosing realistic losses and typical values of xi an d 
k\ (see detailed calculations for a similar cavity configuration in [[Ll]). Assuming xi = 
ki = 3.4 x 10 4 s~ 1 , 71 = 2.4 x 10 8 s _1 , 72 = 3.6 x 10 8 s _1 , we find for the threshold pump 
power P t h = 74mW. For these parameters the power of the cavity output field of u/4 mode 
increases as P° ut = 2ft~7i7ii = 3.4 x 10~ 3 (E/E th — 1)W, while the output power of the u/2 
subharmonic remains constant P 2 ° u ' = 2/2.^72^2 = 3.4 x lO" 3 ^. 

A very important characteristic of multiphoton OPO concerning phase information of 
the generated modes, is displayed in phase-space. As we see from (||), the stationary states 
of u/4 subharmonic have four-fold symmetry. There are four stable states of iO\ mode with 
equal photon numbers, but with four different phases which are: f,-f + f,f+7r,-f + ^f - As 
regards to uo/2 subharmonic mode, it displays two-fold symmetry in phase-space, i.e. there 
exist two states with equal intensities but with different phases *, * + n. 

Let us now show the general character of this symmetry in the frame of a full quantum 
treatment of interaction. It is easy to verify that the Hamiltonian ([!]) as well as the density 
operator p obeying Eq.(^) satisfy the commutation relations 



U = exp \ i— a[ai + ma\a2j ■ 

For the reduced density operators of each of the modes pi( 2 ) = Tr 2 (\)(p), (which are con- 
structed by tracing over one of the modes), relations (0) give 



[H 1 ,U} = [p(t),U] = 



(7) 



with the operator 
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Pi(t),U^)] = [p 2 (t),U 2 (n)} = 0, (8) 

where operators Ui(ip) = exp (iipa\ai^ , i = 1, 2 perform rotation by the angle ip around the 
origin in phase-space of complex stochastic variables ot\ and a 2 , respectively. 

One of the most important conclusions of such symmetries is related to the Wigner 
functions W\ and W 2 of the subharmonic modes ui\ and u 2 which provide a large amount 
of information about the states of the mode and also provide a pictorial view. To see 
this we address to the formula for Wigner function expressed through the density operator. 
Using also formulas (f|) we arrive to the symmetries of Wigner functions written through the 
stochastic mode amplitudes cti and a 2 corresponding the operators a x and a 2 . In the polar 
coordinates r, 6 of the complex phase-space plane X = aH 2 Q * = r cos 9, Y = a ~ a * = rsm9 
these symmetries take the forms: 

W 1 (r,9 + n/2) = W 1 (r,9), W 2 (r,6 + tt) = W 2 {r,9). (9) 

We have obtained this property of Wigner functions using the symmetry of Hamiltonian ([[]) 
and dissipation terms in master equation (^|), assuming that the initial state of the modes 
is a vacuum state. Below, we will present the other confirmation to phase-space symmetry 
in quantum trajectories. 

III. THREE-PHOTON OPO 

The dynamics of subharmonic modes u)\ = u/3 and u 2 = 2u/3 for this system is governed 
by master equation (fj) involving Hamiltonian (0) instead of Hamiltonian ([I]) . Accordingly, 
the classical equations of motion for the complex amplitudes are 

— 71CH1 + X2a^E + 2k 2 a 2 al, 



dt 

72«2 + X2<y*\E - k 2 a\. (10) 



dt 

We list below the semiclassical results for stable values of photon numbers and phases of 
the modes. In the above-threshold generation regime, at E > E t h, where 



is the threshold value of E, the results take the form: 

7i72 / . q r~ 2 ^ 2 



where e — E/E-, 



thi 



m 

,^4,f g+3 ^ZTV, (12) 



$ 2tt , 2$ 2tt 

- + — n , 2 = — -— n, n = 0,l,2. (13) 



In the region e < the stability condition is fulfilled only for a zero amplitude steady- 
state solution a.i = «2 = 0, i.e. the subharmonic field excitation exists at the spontaneous 
noise level. The photon numbers are plotted against the pump amplitude in Fig.2. They 
show bistable hysteresis-cycle behavior in a small domain 1 < e < The dashed parts of 
the curves correspond to unstable solutions. As we see, with increasing E/E t h the photon 
number of the cu/3 subharmonic increases as n\ ~ 8 g^p £ 2 , while the photon number of the 
2c<j/3 subharmonic goes to the saturated value n 2 = 7 2 /9/c|. 

The very important peculiarity of the system proposed is that the threshold value Eth 
depends on the coupling constant \2 which is related to the second-order x^ susceptibil- 
ity in contrary to the scheme cited in |J, where the threshold of the driving field ampli- 
tude is determined by x susceptibility. However, the photon numbers of the generated 
modes depend on both nonlinear coupling constants X2 and k%. For the realistic parameters: 
k 2 = 2.4 x 10 3 s _1 , X2 = 2.3 x 10 3 s _1 , 71 = 0.9 x lO 8 * -1 , 72 = 1.2 x 10V 1 parametric 
oscillations occur above the threshold pump power P t h = 3W and the threshold values of 
cavity-output field powers are P{ h ~ Am W, P 2 th ~ AmW . Eqs. (p^2|) , (|T3|) imply that in the 
above-bistability regime there exist three states for each of the modes which have equal 
intensities n\ and n 2 but different phases: §> § + ^fj § + 4f and ^ — j, ^ — 
respectively. Surprisingly, both the modes uoi = u/3 and uo 2 = 2u/3 have the identical 
threefold symmetry which also follows from the specific three-photon form of Hamiltonian 
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(0). In the full quantum treatment of cascaded three-photon down-conversion this symme- 
try, in particular, is displayed in the Wigner functions. Following the method presented 
above, we have obtained: 

Wi(r,9 + 2*/3) = Wi(r,9), i = l,2. (14) 



IV. QUANTUM SIMULATION: MULTISTABILITIES AND THRESHOLD 

BEHAVIOR 

The second part of this paper is devoted to the study of quantum-statistical properties 
of cascaded multiphoton OPO in a dissipative cavity. Our aim is to analyze the phase-space 
properties of the subharmonic modes in a quantum regime in the presence of quantum noise 
on the basis of the Wigner functions. 

Note that an exact analysis of quantum optical nonlinear systems including dissipation 
effects is a difficult problem which was solved for a few single models (see [ff2l and references 
therein). In this field of research, analytical expression for the Wigner function of two- 
photon OPO was derived in [O. However the quantum-statistical properties of cascaded 
OPO's can not be analyzed using single analytical formulas. Our approach is based on the 



quantum-jump simulation method also known as the state- vector Monte-Carlo method |L4 
This method considers not the density matrix but deals with the state vector | ty a (t)) which 
is a member of the ensemble of state vectors and index a indicates the realizations. The 
procedure adopted in the quantum-jump simulation for the system Hamiltonians ([I]) and (0) 
consists in the following. The evolution of a member of the ensemble of pure states | \I/ Q (t)) 
over a short time 5t is governed by the non-Hermitian effective Hamiltonian if there is no 
quantum jump 

H i(2)eff = #1(2) - ift(7ia{ai + 12^2) ■ (15) 

Such evolution must be completed by the possibility of quantum jumps which change 
the state vectors of both modes. On the whole the state transforms in the following way 
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V a (t + St)) 



(1 - {5tH l{2)eff ) | # a (t))/(l - dp) 1 / 2 , with probability 1 - Sp, 
V^ai | ^ Q (t))/(5j5j/5t) 1/2 , with probability 5p h i = 1,2, 

(16) 



where <5p = <5pi + <5p 2 . 



We generalize the method presented in [15] for the quantum-jump simulation of the 



Wigner functions. Details of analogous calculations for the second harmonic generation can 



be found in |16|.The numerical simulations are performed in the truncated Fock basis of 
both subharmonic modes ui and u 2 in the regime of strong nonlinear coupling between the 
subharmonic modes, i.e. xi,2 ^ 7i,72 and k 12 ^ 7i,72- We note that the OPO with such 
extremely large nonlinearities are not realized in practice and in this part of the Letter we 
do not intend to give results close to an experimental situation but discuss the fundamental 
problems of multiphoton entangled states. 

To be short, we shall give the results for three-photon OPO, where one kind of phase- 
space symmetries is realized. The steady-state Wigner functions for u)\ = u/3 and u 2 = 2u/3 
modes are calculated in the basis of n = 45 photon number states for each mode. The results 
may be written as depending on the dimensionless parameters: e = E/E th , k 2 / r fi, 72/71- 
We calculate the Wigner functions of U\ = uj/3 subharmonic averaged over 1000 trajectories 
in the complex phase-space plane X, Y. It is obvious that below the threshold the Wigner 
function is single-humped and centered at X = Y = 0. This state indicates appearance 
of three arms as the consequence of small entanglement between the modes. Increasing 
the driving field E we enter the bistability domain and observe the occurrence of three 
additional side-humps. They correspond to the above-threshold steady-states (12), (13) 
with equal photon numbers and threefold symmetric phases. We see explicitly that in this 
domain the behavior of our system is actually four-stable. With further increase of E/E t h 
the central hump disappears, while the side-humps increase and we turn to a manifestly 
above threshold oscillation regime when we observe the phase space three-stability. This 
case is shown in Fig. 3 for the parameters k 2 /ji = 0.2, 72/71 = 0.4, e = 1.59. Analogous 
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behavior has the Wigner function of 2c<;/3-subharmonic.. The difference between the results 
of two modes is displayed in the form of squeezing the humps. 

V. CONCLUSION 

Thus we have shown that the idea of using the cascaded photon splitting processes in a 
cavity opens a novel route for the realization of a multiphoton parametric oscillator. Two 
devices are proposed which transform the pump coherent field into two subharmonics at 
frequencies uo/2 and uo/A (configuration A) and uj/3 and 2cj/3 (configuration B). Further- 
more, these models of OPO are to be subject to an experiment. We have been convinced 
that the pump threshold of each of the OPO is only expressed through the second-order 
susceptibility. As a result, comparatively low values of thresholds have been obtained for 
typical parameters. The important features of multiphoton OPO obtained are phase-space 
symmetries and multistabilities of subharmonics. In addition, it is shown that subharmonics 
"remember" the phase of the pump field. Our numerical quantum Monte-Carlo simulations 
display phase-space multistabilities and threshold behavior of subharmonics in the presence 
of dissipation. 
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Figure Captions 

Fig.l. Principal scheme of the cascaded parametric oscillator. The cavity is resonant 
at two frequencies U\ = u/4 and u 2 = oj/2 of the subharmonics. The phase matching 
condition for the process cj^c<j/2 + c<j/2is satisfied in the first medium, and for the process 
ui/2 u/4: + u/4 in the second medium. 

Fig. 2. The normalized mean photon numbers ni/nf 1 (curvel) and r^/n^ 1 (curve 2) of 
uj/3— and 2cj/3— modes versus the scaled pump intensity parameter e 2 = E 2 /E 2 h . The 
threshold values of the photon numbers are n\ h = and nf 1 = j^i- The dashed parts of 
the curves 1,2, and zero-amplitude (curve 3) solutions for both modes describe the unstable 
steady-state solutions. 

Fig. 3. The steady-state Wigner function ofu/3 mode which demonstrate the phase-space 
multistability in the above-threshold range. 
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